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ABSTRACT 

In this paper we pursue the suggestion by Shu, Lubow & Anderson (1979) and Wang (1995) that 
contact discontinuity (DSC) may exist in the secondary in the expansion TRO (thermal relaxation 
oscillation) state. It is demonstrated that there is a mass exchange instability in some range of mass 
ratio for the two components. We show that the assumption of constant volume of the secondary 
should be relaxed in DSC model. For all mass ratio the secondary alway satisfies the condition that no 
mass flow returns to the primary through the inner Lagrangian point. The secondary will expand in 
order to equilibrate the interaction between the common convective envelope and the secondary. The 
contact discontinuity in contact binary undergoing thermal relaxation does not violate the second law of 
thermodynamics. The maintaining condition of contact discontinuity is derived in the time-dependent 
model. It is desired to improve the TRO model with the advanced contact discontinuity layer in future 
detailed calculations. 
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1. INTRODUCTION 

W Ursae Majoris (W UMa) binary stars were first 
' thought to be a particular binary population due to their 
abnormal mass-radius relationship, namely, the so-called 
Kuiper's paradox, R2IR1 = {A'h/Mif-'^^ (Kuiper 1941). 
, These particular binaries appear to consist of two main- 
■ sequence stars that possess photospheres exhibiting the 
almost same effective temperatures for the two compo- 
nents despite the fact that typical mass ratio in a system 
' is 0.5. It was originally proposed that a common convec- 
I five envelope may be formed due to dynamic equilibrium 
(Osaki 1965), and mass and energy transfers would take 
place in CCE in order to interpret the Kuiper's paradox 
(Lucy 1968) although the specific mechanism of energy 
transfer for the circulation has not been fully understood 
' (Robertson 1980, Sinjab, Robertson & Smith 1990, Tas- 
soul 1992). It is now firmly believed that W UMa stars 
are contact binaries in which both components are full of 
their Roche lobe, showing strong interactions (Mochnacki 
1981). Lucy's iso-entropy model (1968) as a zero-order 
model with thermal equilibrium, however, cannot explain 
the color - period diagram by Eggen (1967) which leads 
to the establishments of two parallel first-order theories of 
thermal relaxation oscillation (TRO) and contact discon- 
tinuity (DSC). TRO model was advanced by Lucy (1976), 
Flannery (1976), and Robertson & Eggleton (1977), who 
suggest that contact system can not reach thermal equi- 
librium at dynamical equilibrium configuration and may 
thus undergo thermal relaxation oscillation. DSC theory 
was proposed by Shu and his collaborators [Shu, Lubow 
& Anderson 1976, 1979, Lubow & Shu 1977, also see Bier- 
mann & Thomas (1972) and Vilhu (1973) for some earlier 
elements of DSC model] who hypothesize contact binary 
can attain thermal and dynamical equilibrium but there 
is a temperature inversion layer in the secondary. With 
great attempts the so-called Kuiper's paradox and period 



- colour diagram may be resolved by the two different 
hypotheses independently. However, both theories have 
some difficulties to explain observations, such as the so- 
called W-phenomenon, i.e., the less star is hotter than the 
massive component (Binnendijik 1970) (see a concise sum- 
mary of observations and theory by Smith 1984). Espe- 
cially there are great debases between the two in nature in 
their simplest version (Kahler 1989). Observational stud- 
ies continue and the theoretical controversy still remains 
(Rucihski 1997). These imply that the first-order theories 
of contact binary (i.e. TRO and DSC) should be improved. 

The intensive disputes by the two contending schools 
(Lucy & Wilson 1979, Shu, Lubow & Anderson 1979) lead 
to an intriguing suggestion by Shu, Lubow & Anderson 
(1979), Shu (1980) (from theoretical viewpoint) and Wang 
(1995) (from the analysis of observational data) that TRO 
theory needs the contact discontinuity in some phases. 
Some observations seem to support TRO theory (Lucy 
& Wilson 1979, Hilditch et al 1989, Samec et al 1998). 
Although some criticisms about DSC model exist (Shu, 
Lubow & Anderson 1980), this theory is still attractive 
because it is successful in many aspects (Smith 1984). The 
contact discontinuity may be ironed out within the thermal 
timescale (Webbink 1977, Hazlehurst & Resfdal 1978, Pa- 
paloizou & Pringle 1979, Smith, Robertson & Smith 1980) 
in a steady state, however, the existence of time-dependent 
contact discontinuity can not be excluded because it does 
not violate the second law of thermodynamics (Papaloizou 
& Pringle 1979). However a detailed analysis is needed for 
this. It is highly desired to reconcile the two theories not 
only for removing the discrepancies but also for explaining 
more detailed observations (Shu 1980). 

The difficulties of pure TRO and DSC models in ex- 
planation of W-phcnomenon motivate us to explore the 
possibility to develop a second-order theory. The inter- 
action between the secondary and the common convec- 
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tive envelope is thought as an important role in the W- 
phenomenon, Wang (1994) find that the W-phenonienon 
can be explained by the released gravitational energy of 
the secondary through its contraction corresponding to the 
TRO contracting phase in W-type contact systems. This is 
encouraging and leads to the suggestions by Wang (1995) 
from a sample with 32 contact systems that the A-type 
systems may undergo thermal relaxation oscillation with 
contact discontinuity whereas the contraction of secondary 
in W-type systems irons out this contact discontinuity. 

The over-riding virtue of a contact discontinuity is that 
it gives a clear mechanism for making the secondary physi- 
cally larger, and the primary physically smaller, than their 
main-sequence single-star counterparts, as is needed to sat- 
isfy the Roche-lobe filling requirements of the Kuiper para- 
dox (Lubow & Shu 1977). On the other hand, if the sys- 
tem cannot be maintained in steady state by heat-carrying 
flows, then the capping of the radiative heat flow from the 
secondary by the hotter overlying common envelope should 
lead to an expansion of the secondary, with a resulting 
transfer of mass from the secondary to the primary. Such 
Roche-lobe overflow, from a less massive star to a more 
massive star to a more massive one, is known to occur 
slowly, so the ultimate breaking of contact caused by the 
expansion of the binary orbit takes a relatively long time. 
Once contact has been broken, however, the common en- 
velope disappears; the secondary is no longer capped; and 
it will begin to shrink toward its normal single-star size. 
Conversely, because the larger area of the common enve- 
lope is no longer available to carry away much of the pri- 
mary's interior luminosity, the primary can no longer be 
maintained at its suppressed contact size, and it will begin 
to overflow its Roche lobe. The transfer of mass from a 
more massive star to a less massive one is known to be 
unstable (e.g. Paczynski 1971), and the rapid shrinkage 
of the binary orbit causes the system to come into con- 
tact again. The re-establishment of the common envelope 
and the capping of the secondary results in its refilling its 
Roche lobe. Thus, would the DSC hypothesis provide the 
physical mechanism for the TRO hypothesis, together with 
a justification why the duty cycle is long for the contact 
phase and short for the semidetached phase, as is required 
by the observational statistics. The rest of this paper at- 
tempts to establish the above ideas on a more quantitative 
basis. 

This paper is organized as following: the instability of 
Roche lobe and its operation in contact system are found 
in Sec. 2; the surviving condition of DSC layer is derived 
from the thermodynamics in Sec. 3; and the conclusions 
are remarked in last sectionr. 

2. THE ROCHE LOBE INSTABILITY AND DSC MODEL 

It is generally believed that the two components of W 

UMa stars share an optically thick common convective 
envelope due to the dynamics equilibrimn (Osaki 1965, 
Mochnacki 1981). The redistribution of the total lumi- 
nosities (the plus of luminosities of each star), which takes 
place in CCE, deals with comprehensive fluid processes 
(Lucy 1968). Why and how to redistribute the luminos- 
ties is the main task to theoreticians. The debate of the 
exsiting theories of contact binaries had been attracted 
much consideration between 1970s and 1980s (Lucy & Wil- 
son 1979, Shu 1980, Shu, Lubow & Anderson 1981). Shu 



(1980) clearly stated that the two superficially distinct the- 
ories are complementary with the crucial theoretical issue 
to be resolved being the secular stability of temperature 
inversion layer from his thought-provoking analysis. Here 
we argue that DSC layer is a natural results of TRO theory 
via the mechanism of Roche lobe instability, showing the 
presence of DSC layer during the expansion TRO phase. 

In the following discussions we assume that the total 
mass and angular momentum are conserved, neglecting 
the spin angular momentum of two components. These 
assumptions are basic and the same in TRO theory, but 
they are unnecessary in the DSC theory. In principle, the 
two assumptions put more strong constraints on the the- 
oretical model. In the conserved systems there are mainly 
two other parameters: mass ratio q, and mass ratio chang- 
ing rate due to mass exchange q, to determine the struc- 
ture of the contact binaries in TRO theory. The most 
serious shortcoming of TRO (mentioned in the previous 
section) is a strong indicator that we should relax some of 
assumptions in TRO model. One possible way to remove 
this shortcoming is to supplement the interaction between 
CCE and the component. This inclusion may reconcile the 
two contending schools each other (Wang 1995). 

We first show that the instability of mass exchange may 
prevent from the mass in the secondary being pushed into 
the primary through the inner Lagrangian point Li due 
to the lid effects of CCE placed on the secondary (Shu, 
Lubow & Anderson 1976). For a contact system with to- 
tal angular momentum ^7 in a circular orbit and total mass 
j\4 = Ml -|- M2, the separation between components reads 
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where q = M2/M1 (for the convenience we take g < 1), 
and G is the gravitation constant. The Roche lobe ra- 
dius i?L of the secondary approximates for all mass ratio 
(Eggleton 1983) 
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The Roche lobe of the primary will be obtained when we 
replace q hy 1/q. It is important to note that the Roche 

lobe is changing due to the mass transfer between the two 
components. The variation rate of the Roche lobe due to 
mass transfer between the two components reads 
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and then we have the timescale for this change with the 
helps of dlni^LM = {d In Ri,/dq){dq/dt) 
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where tyi is the timescale of mass transfer defined as 
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here the parameter Mi is the rate of mass transfer, and 
the function f{q) is 



defined as p oc A'h/Rl- Wc thus have 
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The function f{q) represents the ratio of the two 
timcscalcs. Wc have calculated the function f{q) in Figure 
1, showing its value for the range of q from 0.0 to 1.0. If 
f{q) > then Roche lobe will expand with increases of q 
or shrink with decreases of q. If f{q) < then Roche lobe 
will shrink with the increases of q or expand with decreases 
of q. It is very important to address that if |/((?)| < 1 then 
the expansion or shrinkage will be rapid than the process 
of mass transfer in the contact system known from equa- 
tion (4). From Figure 1 the Roche lobe of the secondary 
expands with the increases of mass ratio. The expansion 
timescalc is shorter than that of gaining mass from the pri- 
mary until the mass ratio g > 0.8, indicating that the sec- 
ondary is capable of swallowing more mass when q < 0.8. 
In this range of mass ratio the mass gaining is unstable. 
The Roche lobe of the primary will shrink due to the mass 
exchange. There also is an instability of mass exchange in 
the range of g < 0.35. This implies that the timescale of 
Roche lobe shrinkage due to mass transfer is shorter than 
that of mass transfer. This means that Roche lobe shrinks 
more rapid than the mass loss. This mass exchange in- 
stability plays important role in the structure of contact 
binaries. The maximum of mass ratio for the instability 
of the primary is 0.35. It is believed that the mass trans- 
fer will be more efficient when q < 0.35 with the presence 
of the instability of mass exchange. Thus it is expected 
that the mean mass ratio of A-type systems will be less 
than that of W-types. This is consistent with the observa- 
tions. It should be noted that here wc do not specify the 
mechanism for the energy and mass transfer. Of course 
the direction of mass exchange between the two compo- 
nents is determined by the relative potential of the star 
surface. Here we are trying to show the instability of mass 
exchange, namely, represented by !/(<;) | < 1. 

In the original DSC version the rising convective ele- 
ments interior to the Roche lobe of the secondary cannot 
penetrate into the common convective envelope because 
the resulting buoyancy deficit opposes such penetration. 
SLA76 argued that there is a mass flow pushed by the 
slight excess of pressure due to a slight heating of the in- 
terior of the secondary under constant volume. It is very 
important to note that this mass flow process is based 
on the assumption with constant volume of the secondary. 
The mass flow was estimated by SLA79, however, their 
estimation follows up another assumption that all the en- 
ergy transferred from the primary to the secondary radi- 
ates again from the secondary, neglecting the interaction 
between CCE and the secondary. Now we can work out 
a condition inhibits the returning of mass flow to the pri- 
mary through the inner Lagrangian point Li. If there 
is no excess pressure, then the mass flow stops. This is 
equivalent to dp/dt < if the contact discontinuity be- 
ing lower than the temperature of CCE survives, namely, 
T =constant (we assume the gas is ideal), where p is the 
mean mass density within the Roche lobe of the secondary 
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then the condition no returning of mass flow reads 



(7) 



(8) 



We draw the line = Sqtu in Figure 1. It is obvi- 
ous that all the value of f{q) is always less than 3q. This 
means that all the cases satisfy the condition that no mass 
flow returns to the primary even beyond the mass ex- 
change instability. Therefore the assumption of constant 
volume of the secondary should be relaxed in the advanced 
DSC model. The contact discontinuity is time-dependent 
from this viewpoint at least, coinciding with that DSC 
layer could be maintained in a time-dependent model (Pa- 
paloizou & Pringle 1979). 

3. THERMODYNAMICS OF DSC LAYER 

By defining the thermal timescale as ixh = 
J^_g^{4Trr'^ pvc)~^dm, Webbink(1977) first showed that 
the thermal diffusion time scale in the common envelope 
is typically of the same order as the dynamical timescale 
(is roughly of one orbital period). This makes the contact 
discontimuty disappear within one orbital period. Wc call 
the thermal diffusive process as interaction e. Papaloizou 
& Pringle (1979) show the steady contact discontinuity 
violates the second law of thermodynamics, but the time- 
dependent contact discontinuity may exist. However in the 
time-dependent model it is the interaction e that keeps the 
contact discontinuity in contact system undergoing ther- 
mal relaxation oscillation. The controversy of inner struc- 
ture may be removed by this kind of interaction (Wang 
1995). With the helps of the conservation of mass and 
momentum we can rewrite the energy equation beyond 
the energy generation region as 



p|(* + T.) 



-pTv-Vs-V -F + e, 



(9) 



for the inviscid fluid (e.g. Webbink 1977, and 1992 ApJ, 

396, p378 for the crratmn), where t is time; p, the den- 
sity; T, the temperature; v, the velocity; s, the specific 
entropy; F, the energy flux radiated from the star; e, the 
energy density absorbed by the secondary in the unit time 
due to the interaction with CCE; and the gravitational 
energy per unit mass. Following the assumption by Shu, 
Lubow & Anderson (1980) that the specific entropy s can 
be decomposed in terms of a barotropic and a baroclinic 
one as s = so(5'd) + si{x, t), wc integrate the above equa- 
tion over the volume enclosed by the equipotential surfaces 
C and D, and obtain 
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where S = § psdV, AM = J pdV, AV = g dV is the 
volume enclosed by the two surfaces C and D, dA is the 
area of the surface of contact discontinuity, and n is its 
normal vector. For the time-dependent case we assume 
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that the last two terms offset approximately as in steady 
case (Shu, Lubow & Anderson 1980), thus we have more 
physically concise form of equation (10) 

The first term of right hand in equation (11) represents 
the entropy increases due to mass exchange between CCE 
and the secondary, the last term does the same meanings 
but due to energy exchange, the second term does the 
entropy decreases of entropy due to the Roche lobe expan- 
sion. The enclosed volume is an open system undergoing 
mass and energy exchanges with its surroundings rather 
than an isolated volume. This equation also tells us the 
resulting expansion due to the interaction e if the contact 
discontinuity layer survives: 1) if there is only exchanges 
of energy by thermal diffusion, namely, AM =const, we 
have 

This clearly states that the surviving of contact disconti- 
nuity must be provided by the expansion. Detailed calcu- 
lation should be done in the future. 2) There is a mass 
exchange between CCE and the secondary accompanying 
the energy interaction, i.e., d{AM)/dt > 0, the expansion 
is at least 

This equation predicts the secular change of orbital period 
due to the shift of mass ratio. 3) If the secondary keeps 
constant volume as originally suggested by Shu, Lubow & 
Anderson (1976), the term d'i'/dt = 0, we always have 
dS/dt > which means the discontinuity will be ironed 
out with the thermal diffusive timescale. The only possible 
way to relax this condition is the inclusion of the changes 
of Roche lobe. This way will permits us unifying the two 
contending hypotheses. According the simplest version of 
star structure, equations (12) and (13) will provide the 
expansion velocity 

.«>«.(y(f)+r.,(*„te.!M, (14) 

where e = (GM2AM/R2)/V2 is the mean density of the 
gravitational energy between the secondary and CCE, 
.92 = GAl2 1 -R| the gravitational acceleration, V2 is the vol- 
ume of the secondary. One should note that in the above 
estimation we neglect the down directed propagation of 
energy to the interior with dynamical timescale. There- 
fore the present estimation is somewhat higher than that 
of the actual. Both of the mass exchange and energy in- 
teraction result in the expansion of the secondary, we thus 
have the minimum velocity 

y = max(wint,URL), (15) 

where = dR\,jdt is the expanding velocity of the 
Roche lobe. This is the condition maintaining contact dis- 
continuity. Prom the viewpoint of total energy (by the 
nuclear) conservation, the exhausted energy (i.e. e) to ex- 
pand the secondary lowers the re-radiating efficiency of 
the transferred energy from the primary. The lower tem- 
perature of the secondary than the primary may be an 
indicator of the presence of contact discontinuity. 



4. CONCLUSIONS AND DISCUSSIONS 

Introducing a discontinuity of temperature by Shu and 
his collaborators (1976, 1979) the thermal instability of 
binary (Lucy 1976, Flannery 1976) can be suppressed, but 
its maintenance of DSC layer opens. In this paper we 
try to construct the physical scenario of time-dependent 
model of contact binary. Only two assumptions that total 
mass and angular momentum of the contact system are 
conserved are employed in this paper. It is found that the 
mass exchange results in the instability of Roche lobe in 
some ranges of mass ratio. We show that this instability 
always satisfies the condition that keeps the mean density 
of the secondary dlnp/dt < 0. Therefore it ensures that 
no mass flow returns to the primary through the inner La- 
grangain point Li . The second-order theory predicts that 
the contact binary may be in oscillations take place about 
a state with a contact discontinuity. The temperature dif- 
ferences of DSC layer across the interface is determined by 
the expansion velocity. 

The existing TRO and DSC theories (Lucy 1976, Flan- 
nery 1976, Robertson & Eggleton 1977; Shu, Lubow & 
Anderson 1976, 1979) neglect the effects of interaction e 
between CCE and the secondary. Here we argue that it 
is the contact discontinuity layer that results in the inter- 
action between CCE and the secondary and in the mean- 
while it is the interaction that maintains the contact dis- 
continuity. We find this interaction e can result in some 
interesting issues. First it is the reason why the tempera- 
ture of the secondary in A-types is lower than that of the 
primary. Second the maintenance of contact discontinu- 
ity needs faster mass transfer which breaks down the deep 
contact. Thus the shortcoming of TRO theory will be re- 
moved. It is highly desired that the unification of TRO 
theory and DSC model should be calculated in order to 
discover the nature of the contact binaries. 

In the present work we do not specify the mechanism 
of thermal diffusion process. Although we have not per- 
formed the time-dependent model unifying TRO and DSC 
hypotheses, this time-dependent unified theory might give 
some predictions. First the secular behavior of period 
change of A-type systems is violent than that in W-type 
system in order to survive the existence of contact discon- 
tinuity. Second, the maintenance of contact discontinuity 
may lead to the radial oscillation of the secondary with 
period from a few to several ten minus. The interaction 
between CCE and the secondary drives such a oscillation 
similar to the k- mechanism working in other types of stars. 
It is thus expected to find the light variation during the 
primary eclipse as another probe of contact discontinuity 
in A-type systems. 
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